Proposal for a ferromagnetic multiwell spin oscillator 
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The highly nonlinear coupling of transport and magnetic properties in a multiwell heterostructure, 
which comprises ferromagnetic quantum wells made of diluted magnetic semiconductors, is theo- 
retically investigated. The interplay of resonant tunneling and carrier-mediated ferromagnetism in 
the magnetic wells induces very robust, self-sustained current and magnetization oscillations. Over 
a large window of steady bias voltages the spin polarization of the collector current is oscillating 
between positive and negative values, realizing a spin oscillator device. 



Magnetic resonant tunneling structures are prominent 
spintronic devices |l|, which are proposed for spin valves, 
spin filtering 0-Ia], or for realizing digital magnetoresis- 
tance [7|. In nonmagnetic multiwell heterostructures in- 
teresting dynamic effects such as the formation of electric 
field domains and the motion of charge dipoles through 
the structure can be observed [3, Q- Recently, it has 
been predicted that in structures with paramagnetic wells 
these phenomena can be controlled by an external mag- 
netic field |1CHL2J . The paramagnetic systems can behave 
as a spin oscillator device, in which the moving charge 
dipoles generate spin-polarized current oscillations in the 
MHz range [12,-14] . In ferromagnetic multiwell structures 
made of dilute magnetic semiconductors (DMS) in addi- 
tion also the magnetic order will become dynamic |15| . 
This additional degree of freedom enriches the dynamic 
complexity, which up to now is still largely unexplored. 
DMSs |16l - [l8| are made magnetic by doping them with 
transition metal elements, e.g., by incorporating Mn in 
a GaAs crystal host. In several experiments ferromag- 
netism has been generated in the bulk by tailo ring the 
actual particle density electrically or optically [1^, |20| . 
In 2d-confined systems made of DMSs the magnetic or- 
der depends strongly on the local spin density, which is 
influenced by the tunneling current |2l| - [23| . This inter- 
play leads to a novel mechanism of generating robust 
self-sustained current oscillations. 

In this article a detailed study of the dynamics in a 
ferromagnetic four- well system, as a prototype of the ex- 
pected multiwell functionality, is provided. The carrier 
dynamics is described by a self-consistent sequential tun- 
neling model which includes momentum and spin relax- 
ation inside the wells. The feedback effects of both the 
carriers Coulomb interaction and the magnetic exchange 
coupling with the magnetic ions are described within a 
mean-field picture. We give an qualitative explanation 
for the occurrence of spin-polarized current oscillations, 
which realizes the functionality of a spin oscillator device. 

The band profile of a four-well heterostructure with 
two ferromagnetic quantum wells made of a DMS, e.g., of 
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FIG. 1: (Color online) Schematic scheme of the band profile 
of the ferromagnetic four-well structure. The exchange inter- 
action of the magnetic ions in the second and third well is 
mediated by the carriers tunneling in and out of the wells. 



GaMnAs, is sketched in Fig.[TJ In the ferromagnetic wells 
the ferromagnetic order of the magnetic ions is mediated 
by the itinerant carriers. Describing the mutual exchange 
interaction within a mean-field picture allows to derive an 
analytic expression for the exchange splitting A° of the 
subbands in the ferromagnetic wells [l|, [2l|, [22, 124 
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where Jpd is the coupling strength between the impurity 
spin and the carrier spin density (in case of GaMnAs p- 
like holes couple to the d-like impurity electrons), z is the 
longitudinal (growth) direction of the structure, 7iiinp(-z) 
is the impurity density profile of magnetically active ions, 
-00(2:) labels the well ground wave function, and s = 1/2 
is the particles spin. The Brillouin function of order S 
is denoted by Bs, where 5' is the impurity spin, which 
for Mn equals 5/2. We assume a homogenous impurity 
distribution making the spin polarization ^i — s{nl — nj) 
the determining factor for the exchange splitting A^. 

The well spin polarization can be rapidly changed by 
in- and out-tunneling particles. However, the magnetic 



impurities need some time to respond until the corre- 
sponding mean field value A° is established. In the case 
of GaMnAs, experimental studies of the magnetization 
dynamics revealed typical response times ta of about 
100 ps [25 ■ We model the magnetization evolution in 
the ferromagnetic wells i = (2, 3) within the relaxation 
time approximation, dA^/di = — (A^ — A^)/rA. 

The transport through weakly coupled multiwell struc- 
tures is well described by a sequential tunneling model, 
which is presented in detail in Refs. P,'l2!|. The transport 
model includes the following assumptions: (i) The wells 
are weakly coupled and tunneling can be described by 
an transfer Hamiltonian formalism, (ii) The particles in 
each well always stay in a quasiequilibrium state due to 
elastic and inelastic scattering, which is assumed to con- 
stitute the shortest time scale in the problem (rscatt ~ 1 
ps). (iii) Particles in the first excited subband relax 
rapidly to the ground level by inelastic scattering pro- 
cesses, e.g., due to emission of LO-phonons. This hap- 
pens within a few ps, which is much smaller than the 
other relevant time scales occurring in the problem, (iv) 
The spectral function of the subbands is a Lorentzian 
with the broadening given by 7 = /i/2Tscatt- (v) Dur- 
ing tunneling processes the parallel momentum and spin 
are conserved, i.e., effects of interface roughness are not 
taken into account, (vi) Spin flipping processes inside 
the well can be described by a single spin relaxation time 
Ts ■ Based on these approximations rate equations for the 
spin-dependent particle densities nf , [a =t,i= =tl/2) in 
the ith well can be found 



d^r _ ja _ JO 

At ^^-1-' '^^ 



0,cr 



>j+l 



(2) 



where the tunneling currents Ji^i^i are derived mi- 
croscopically within the Kubo- formalism [3, l26|. The 



0,(T 



quasiequilibrium particle spin densities n^ ' are calcu 



a 



,o,t 



^0,^ 



o,t 



0,4, 



nl' ' — n^'" = ^j^/s, with ^? resulting from inverting Eq. [T] 
for the actual A^. 

Finally, the electron-electron interaction is taken into 
account within a discretized Poisson equation, yielding 
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with Fi denoting the electric field at the ith barrier, e is 
the elementary charge, e denotes the dielectric constant, 
and TVback labels the background charge. 

The tunneling currents are strongly sensitive to the rel- 
ative alignment of adjacent subbands and since the quan- 
tum well levels Efj = Ej + (j)i — aAi explicitly depend 
on the magnetic exchange splitting and the electrostatic 
potential (pi all equations are coupled highly nonlinearly. 

For the numerical simulations we used generic parame- 
ters representing GaAs and GaMnAs quantum wells, re- 
spectively, m — 0.5 mo, Sr — 12.9, Vbar = 300 meV, 
d = 20 A, w = 20 A, £^0 = 80 meV, Ei = 188 meV, 
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FIG. 2: Current- Voltage (IV)-characteristics of the four well 
structure. Self-sustained current oscillations appear in the 
voltage range of V" = 0.22 V to V' = 0.42 V, in which both 
the maximum and minimum values of the oscillations are in- 
dicated by the solid and dashed lines, respectively. 



fJ- 



nm^ [2||, 7 



100 meV, riimp = 0.5 x lO^Ocm-^, Jp^ = 0.15 eV 
2 meV, and Ts =0.1 ns, where m is the 
effective mass with mo denoting the free electron mass, 
and Sr is the relative permittivity, Vbar is the barrier 
height, d and w are the barrier and quantum well widths, 
Eq and Ei are the energies of the two lowest subbands, 
and fie is the emitter Fermi energy. The background 



charge iVback = 0.1 nj, 



10^^ cm ^ is considered to 



be only of about 10% of the nominal Mn doping den- 
sity [27] and the lattice temperature is T = 4.2 K. The 
boundary currents Jo^i and Jn-^n+i with N denoting 
the total number of quantum wells are determined by 
the tunneling from or into a three-dimensional Fermi sea 
at the emitter and collector side, respectively. For these 
currents analytic expressions can be derived [8|, |26| . 

Figure [2] shows the current-voltage {I — V) characteris- 
tics of the investigated structure. We obtain the typical 
sawtooth-like structure, as can be expected for a multi- 
well heterostructure. However, in a broad voltage range 
from V = 0.22 V to I^ = 0.42 V the current does not 
reach a steady state value; instead self-sustained current 
oscillations occur. In this region also the spin polariza- 
tion of the collector current Jn~¥N+i is oscillating be- 
tween positive and negative values realizing in this way a 
spin oscillator device. The effect could serve to generate 
ac spin currents for investigating ac spin injection phe- 
nomena. The maximum amplitude is found at T^ = 0.32 
V, as illustrated in Fig. [31 Along with the current, also 
the well magnetization is oscillating, as shown in Fig. 2] 
The phase shift between the magnetizations (one well 
points into up-direction while the other is down) is al- 
most independent of the applied bias. 

By analyzing the local spin densities and the flow- 
ing tunneling currents in-between adjacent wells a qual- 
itative understanding for the underlying mechanism of 
these oscillations can be found: (i) In the bias region, in 
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FIG. 3: Spin polarization of the collector current versus ap- 
plied bias V in the voltage range in which oscillations are 
occurring. The majcimum and minimum values of the oscilla- 
tions are indicated by the solid and dashed lines, respectively. 
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which self-sustained oscillations are occurring, the cur- 
rents Ji_>.2 and J2^3 are governed by tunneling from the 
ground state to the first excited subband of the neigh- 
boring wells, as shown schematically in Fig. [TJ (ii) Ac- 
cording to an initial small positive energy splitting, spin 
up particles accumulate in the second well and, hence, 
A2 increases, (iii) During this accumulation, particles 
start to tunnel to the third well. At some point the spin 
extraction to the third well overwhelms the spin injec- 
tion from the first well. Therefore, with time the spin up 
particles in the second well are completely extracted and 
spin down particles become the majority species therein. 
This spin inversion comes along with an inversion of the 
exchange splitting A2 oc (n^ — nj,) from positive to neg- 
ative values, (iv) Now the whole process starts from the 
beginning with the role of spin up and spin down particles 
being interchanged. 

The frequency of the oscillations is of the order of 
about 0.25 GHz and depends strongly, as can be ex- 
pected, on the particle tunneling time rj, the spin relax- 
ation time Tg, and on the response time of the Mn-ions 
ta . Since ta is usually the largest time-scale in the prob- 
lem, it determines the final frequency of the oscillations. 
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FIG. 4: (a) Exchange splitting Ai in the second (solid line) 
and third (dashed line) well versus time £ at V = 0.32 V. 
(b) Spin polarization of the collector current versus time t at 
V = 0.32 V. The oscillations frequency is 0.25 GHz. 



In summary, we have shown that in ferromagnetic mul- 
tiwell structures an interesting dynamic interplay of res- 
onant tunneling and carrier mediated ferromagnetism 
emerges, which results in stable self-sustained current 
oscillations. The underlying mechanism based on spin 
inversion is completely different from the one in para- 
magnetic multiwcll structures, in which moving charge 
dipoles are formed. The ferromagnetic heterostructure 
may be useful as a spin oscillator device. 
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